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A  RAPID  COMPUTATIONAL  MODEL  FOR  ESTIMATING  THE 
PERFORMANCE  OF  COMPLIANT  AIRFOILS  IN  CASCADES 

F.  Sisto  and  M.  Avila 
SUMMARY 

By  designing  the  blades  in  a  turbomachine  to  have  a  specific 
schedule  of  structural  stiffness  (typically  more  compliant  than 
normal)  it  is  possible  to  endow  the  end  product  with  an  aero- 
thermodynamic  performance  characteristic  that  is  different  from  the 
rigid-bladed  counterpart.  Presumably  the  resulting  characteristic 
of  the  compliant  turbomachine,  within  certain  operational  limits, 
is  preferable  to  the  more  rigid  turbomachine.  Since  there  is  no 
active  controller  involved  in  changing  the  blade  geometry  in 
response  to  changing  operating  conditions,  the  process  may  be 
thought  of  as  passive  control  by  aeroelastic  tailoring.  In  fact, 
since  the  objective  function  is  a  change  in  the  performance 
characteristic,  the  overarching  process  may  be  termed  "aerothermo- 
elastxc  tailoring":  the  principles  of  static  aeroelasticity  are 
applied  to  a  turbomachine  to  change  its  performance  in  a  thermo¬ 
dynamic  sense. 

The  problem  posed  is  extremely  complex  stemming  from  the  need 
to  couple  the  flow  computation  over  a  deformable  surface  with  the 
structural  dynamics  computation  of  the  body  occupying  the  other 
side  of  the  internal  boundary.  If  no  aeroelastic  instabilities  are 
encountered  at  a  particular  operating  point  the  resulting  con¬ 
figuration  is  fixed  in  time  and  the  static  solution  (flow  and 
structure)  is  a  special  case  of  the  more  general  dynamic 
formulation. 

In  this  situation  it  is  useful  to  develop  an  approximate  model  for 
use  in  the  rapid  exploration  of  many  static  aeroelastic  configura¬ 
tions,  described  by  a  host  of  controlling  parameters.  In  this 
fashion  those  regions  of  parameter  space  showing  potential 
advantages  for  the  compliant-bladed  turbomachine  may  then  be 
computationally  refined  using  more  accurate,  but  more  time- 
consuming,  CFD  coupled  to  discretized  structural  dynamic  codes 
(FEM) .  It  is  the  purpose  of  the  present  model  to  permit  this  rapid 
exploration  of  parameter  space  to  delineate  the  more  promising 
regions  for  the  subsequent  detailed  studies. 

I.  AER.OPYNMIQ,? 

For  the  rapid  computation  of  the  aerodynamic  loads  acting  on 
airfoils  in  cascade,  a  two-dimensional  incompressible  potential 
flow  model  is  chosen  using  a  vortex  sheet  representation  of  the 
thin  airfoils.  Under  these  assumptions  the  upwash  induced  along  the 
chord  of  the  reference  (zeroth)  airfoil  is  given  by 


DTTC  QUAL'T 
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rso  5 


(1) 


fgU)Fa~x)dt 


Here  the  variables  x(and  i)  are  chordwise  coordinates  made  dimen¬ 
sionless  w.r.t.  the  chord,  c,  and  g(x)  is  the  unknown  strength  of 
the  vortex  sheet  representing  each  airfoil.  The  first  integral 
represents  the  contribution  of  the  reference  airfoil  and  the  second 
integral  represents  the  contribution  from  the  remainder  of  the 
cascade  through  the  "influence  function"  F  under  the  assumption 
that  every  blade  has  the  identical  distribution  g(x). 

Figure  1  shows  the  geometric  configuration  under  analysis.  Eq 
(1)  is  derived  in  Appendix  A,  as  well  as  the  inversion  of  the 
integral  equation  for  the  solution  g(x)  under  appropriate  non¬ 
penetration  boundary  conditions  and  specified  vector  mean  velocity 
of  flow  (V, q) . 

The  solution  is  obtained  under  the  classical  thin  airfoil 
coordinate  transformation 


x-Vid-cosQ) 

(2) 

6-cos'1  (l-2x) 

( 2  1  ) 

followed  by  a  discretization  of  the  transformed  variable 

0-—  ,Oitin  (3) 


where  n  is  the  number  of  segments  along  the  chord,  or  n+1  is  the 
number  of  gridpoints  at  which  the  boundary  conditions  are  satisfied 
and  also  the  values  of  g(6)  are  determined. 

Using  results  first  obtained  in  (Sisto,  1967)  the  distributed 
circulation  can  be  expressed  in  the  following  form 


g  (x)  -v. 
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(thus  rendering  f  dimensionless)  and  the  solution  for  f  in  matrix 
format  appears  as  follows: 


f  (0)  [A^]  -Hup } 


(5) 


2 


Here  &  is  the  (n+l)x(n+l)  aerodynamic  matrix  and  U  is  the  ypwash 
column  matrix  or  vector. 

The  expression  for  the  upvash  U  depends  on  the  mean  flow  angle 
a  and  y*  the  camberline  slope  distribution.  The  expression  for  &  is 
quite  complex.  It  reflects  the  geometry  of  the  cascade:  s/c,  0  and, 
again,  y'  as  a  result  of  the  method  of  inverting  the  fundamental 
integral  equation.  The  relevant  expressions  appear  as  Eqs.  (A2) 
(A14)  and  (A16)  determining  &  and  Eqs.  (A14)  (A19)  (A20)  and  (A21) 
for  SJ.  Both  have  been  programmed  for  solution  in  the  Matlab 
language  and  form  part  of  the  program  testoo.ra  and  solve. m 
appearing  as  Appendix  C. 

II.  STRUCTURE 

Each  of  the  identical  airfoil  sections  comprising  the  cascade 
are  modelled  as  two-dimensional  plates  of  small  curvature  and 
variable  thickness.  Each  plate  may  be  thought  of  as  the  character¬ 
istic  or  representative  section  of  an  actual  blade  at  75  or  80 
percent  of  the  span  of  a  typical  cantilever  blade.  As  such  the 
elastic  displacements  may  be  expressed  in  terms  of  the  bending 
moment,  M  according  to 


-  mu)  (6) 

dx2  Eix) 

where  E  is  the  local  stiffness  related  to  Young's  modulus,  cross- 
sectional  moment-of-inertia  and  Poisson's  ratio  if  a  solid 
orthotropic  material  is  under  consideration.  For  an  aeroelastically 
tailored  structure  it  is  presumed  that  E  is  under  the  control  of 
the  designer  within  limits  determined  by  strength  considerations, 
stability,  etc. 

Recognizing  that  changes  in  stagger  related  to  the  untwist  of 
a  three-dimensional  blade  may  be  an  important  aeroelastic  para¬ 
meter,  the  two-dimensional  counterpart  model  will  be  supported  at 
the  leading-edge  and  near  the  trailing  edge  by  springs  whose 
stiffness  may  be  specified.  In  actual  construction  of  a  compliant 
blade  these  springs  represent  the  restraint  to  bending  provided  by 
the  spanwise  spars  required  to  support  the  more  compliant  material 
comprising  the  airfoil  shaped  profile. 

The  elastic  bending  moment  M(x)  is  therefore  a  function  of  the 
aerodynamic  loading  whereas  the  elastic  displacements  depend  as 
well  upon  the  strength  of  leading  edge  and  rearward  springs.  The 
structural  configuration  is  sketched  in  Fig.  2  where  only  a  single 
airfoil  need  be  shown,  the  others  being  identical.  (The  case  of  a 
shrouded  cascade  of  airfoils  would  need  additional  modification  of 
the  model) . 

The  aeroelastic  coupling  is  provided  by  this  aerodynamic 
loading  acting  or.  the  compliant  structure.  (There  may  be  a  small 
amount  of  camber  in  the  unstrained  state  upon  which  the  elastic 
displacements  are  superposed) .  The  aerodynamic  load  is  provided  by 
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dL  -  f VdT  -  CVg(x)d(cx) 


(7) 


dL  -  CV2^ 


— ~  fix)  d(cx) 
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It  is  this  loading  which  is  applied  to  the  plate;  the  details  are 
developed  in  Appendix  B.  It  may  be  noticed  that  dL  has  a  square 
root  singularity  at  the  leading  edge.  However,  under  the  coordinate 
transformation  given  by  Eq  (2),  the  singularity  is  integrable  in 
the  determination  of  lift  and  moment,  viz 

dL  -  gc(l  +  cos9)  f  (0)  d©  (71) 


where 


q  -  V2CV2  (8) 

is  the  dynamic  presssure  based  on  the  vector  mean  velocity,  V.  For 
known  aerodynamic  loading  and  structural  stiffness,  Eq  (6)  may  be 
solved  for  the  elastic  displacements  Y(x).  Subsequently  the  un¬ 
strained  slopes  may  be  augmented  by  the  elastic  slopes  to  obtain 
the  overall  slope  distribution.  This  total  slope  function  comprises 
the  input  to  the  aerodynamics  subroutine  for  computation  of  the 
aerodynamic  loading. 

III.  COMPUTATIONAL  LOGIC 

Conceptually  the  problem  at  hand  is  an  eigenvalue  problem 
similar  in  concept  to  the  problem  of  wing  divergence.  However, 
owing  to  the  complexity  of  the  aerodynamic  formulation,  the  method 
of  solution  is  by  iteration.  An  assumed  blade  shape  produces  a 
specific  aerodynamic  load.  This  load  applied  to  the  blade  structure 
produces  a  specific  schedule  of  displacements.  These  displacements 
define  a  new  blade  shape  and  completes  the  first  loop  of  the 
iteration.  If  the  iterations  converge  the  stable  static  solution  is 
obtained.  If  the  iterations  do  not  converge  the  condition  of 
"divergence”  states  that  the  initial  set  of  aerodynamic  and 
structural  parameters  are  in  supercritical  region  of  parameter 
space  and  no  stable  solution  is  obtainable. 

The  looping  concept  is  illustrated  in  Fig.  3,  a  flow  chart  of 
the  computational  logic  used  to  iterate  the  solution. 

This  chart,  along  with  the  actual  program  statements  (ex¬ 
pressed  in  MatLab  language)  comprise  Appendix  C. 

IV.  PARAMETRIC  STUDY 

A  preliminary  use  of  the  completed  program  demonstrated  was  in 
the  brief  parametric  study  described  below.  The  program  solve. m 
(with  input  provided  from  a  separate  file  testOO.m)  takes  the 
converged  result  from  the  loop  and  computes  a  number  of  related 
parameters. 
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In  order  to  assess  the  proximity  to  stall  a  very  elementary 
approach  is  taken.  The  airfoil  cascade  is  assumed  to  develop  the 
deflection  (lift)  and  drag  provided  by  Howell's  crude  and  early 
correlation  of  the  cascaded  airfoil.  Thus  the  camber  shape,  stagger 
and  solidity  can  be  used  to  enter  Howell's  charts  and  predict  the 
profile  drag  and  lift  (deflection)  coefficients. 

These  data  are  developed  for  the  representative  section  of 
single  rotor  compressor  at  the  radius,  or  spanwise  location,  of  the 
representative  or  characteristic  section.  For  this  study  a  design 
point  wheel  speed  u=468  ft/sec,  stagger  angle  0=45°  and  unit 
solidity  are  chosen*  and  various  graphs  of  loading  coefficient  f 
and  rotor  efficiency,  n ,  can  be  plotted  against  flow  coefficient, 
0.  These  quantities  derived  from  Howell's  correlation  are  not 
presented  as  being  quantitatively  accurate.  However,  they  do 
represent  a  consistent  system  of  predicting  drag-related 
performance  para-meters,  and  the  relative  values  are  thought  to 
indicate  the 

proximity  to  stall,  or  extent  of  separation.  Thus,  the  rotor 
efficiency  at  the  same  flow  coefficient  is  a  measure  of  stall 
mitigation  if  the  efficiency  of  the  compliant  rotor  is  higher  than 
the  efficiency  of  the  rigid  rotor. 

In  Fig.  4  and  the  charts  following  the  results  of  an 
alternative  study  are  shown  for  the  rigid  and  compliant  rotors 


having  the  following  characteristics 

Rigid  Blades 

Compliant  Blades 

Unstrained  Camber,  deg. 

22 

15 

Unstrained  Stagger,  deg. 

45 

45 

Solidity 

1.0 

1.0 

of  Elasticity 

1 . 5xl07psi 

1 . 0xl06psi 

Design  wheel  speed,  fps 

468 

468 

In  this  study  the  compliant  and  rigid  blades  have  been  tailored  so 
as  to  have  virtually  the  same  performance  at  the  wheel  speed,  u=468 
fps.  There  the  system  is  assumed  to  have  an  important  second 
operating  condition  at  10%  overspeed,  or  U=515  fps.  The  performance 
at  this  second  condition  then  discriminates  between  the  two  rotors 
vis  a  vis  stall  mitigation. 

Since  the  aerodynamics  formulation  is  inviscid  in  the  interests  of 
rapidity  of  computation,  an  alternative  assessment  of  the  propen¬ 
sity  to  stall  is  made  on  the  basis  of  several  measures:  incidence 
and  a  leading  edge  loading  parameter  taken  in  conjunction  with  the 
camber  angle.  The  last  is  expected  to  change  significantly  with 
changing  operating  conditions  for  the  compliant  blade. 
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(lower)  for  the  compliant  rotor  at  all  flow  coefficients  in  the 
overspeed  condition.  It  should  be  recalled  that  both  rigid  and 
compliant  rotors  have  been  configured  to  have  virtually  identical 
geometry  and  performance  at  design  speed. 

The  difference  shown  in  Fig.  4  has  not  been  optimized,  or 
maximized.  By  reducing  the  rigidity  of  the  compliant  blades  still 
further,  and  increasing  their  unstrained,  or  "as  fabricated", 
camber,  the  leading  edge  loading  differences  noted  in  the  figure 
could  be  magnified.  However,  this  optimization  process  has  not  been 
carried  forward  at  the  time  of  completing  this  Final  Report.  Only 
the  potential  of  the  method  has  been  demonstrated. 
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APPENDIX  A 


CASCADE  AERODYNAMICS 

Initially  we  wish  to  derive  the  steady  aerodynamic  load  on  a  cascade 
of  camberlines  in  incompressible  flow.  The  usual  thin  airfoil  assump¬ 
tion  is  made  that  the  vorticities  (representing  the  jump  in  tangential 
velocity  across  the  airfoils)  are  distributed  along  the  chordlines  and 
the  boundary  condition  are  satisfied  along  these  same  straight  segments 


/VL  -  -  Z 


Fig.  1  Cascade  Geometry 

On  the  x-axis  the  velocities  induced  by  the  individual  airfoils  are 
given  by* 


c 


★ 

The  Biot-Savart  Law  expresses  the  conjugate  velocity  u,-iv  at  a  general 
field  point.  Subsequent  applications  of  the  Plemelj  formulae  supply 
the  limiting  values  of  u-iv  as  the  camberline  is  approached  from  above 
(+)  or  below  (-) . 
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where  £  =  ^  •+-  A 

pitch,  and  so  on  for  each  airfoil 


/(irr-s) 

and  A  —  S  e  is  the  comple: 


Summing  the  contribution  of  all  the  (infinite  number  of)  airfoils  re¬ 
sults  in 
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Now  a  well-known  result  is  fw  +  L/nirJ  =•  Cl£r+{v  ^  so  that 
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where 


Dropping  the  prime  notation  this  may  be  put  in  the  standard  form 


(0 


where  pr IT  C*f(i  /?)  c/s’  ^ 

F(f"X}T)  "2  "Pc^Hu  \pC^-y)3  —  .  (2b) 

and  1~  (o)  ~  O  i  (y  (o)  ~  O 
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Since  g  is  purely  real,  separation  of  Eq.  (2)  into  real  and  imaginary 
parts  yields 

(i 

^  =  2V  £  3  -0;  +  i* l%) F  <*-*>  £,5  < 

The  boundary  condition  expresses  the  fact  that  the  flow  direction  must 
be  tangential  to  the  blade  surface  (camberline  in  this  case) : 

^  =,  hT  4-  Va^ 

where  V  and  c<  refer  to  the  mean  flow  and  are  defined  in  the  figure*. 
Introduction  of  Equations  (3)  and  (4)  in  Equation  (5)  results  in  the 
following  integral  equation 


+  =iL|g(5)  dL.  (6) 

Solution  of  Integral  Equation 

If  the  LHS  of  this  equation  is  presumed  defined,  then  the  inversion  is 
well  known,  and  is  given  by 

%  W 1  ‘  |  •<  +  ^  Op  [-  f  (->1  -51  ^'(<)  g  (yt))  ^ 

which  may  be  written 

3  Oh  =  1S(X')+  |j(T)  (t> 


In  the  application  of  the  B.C.  the  term  +  b  g(x)  is  omitted  from  M,(x) 
because  of  a  fundamental  approximation  in  thin  airfoil  theory*. 

n - - —  — “ 

Despite  any  local  nonzero  slope  of  the  camberline  the  component  u.  alone 

changes  from  one  side  of  the  boundary  to  the  other.  However,  this  jump 
in  U.  is  the  local  magnitude  of  the  vorticity  and  physically  the  vorticity 
represents  the  jump  in  that  component  of  the  velocity  which  is  parallel 
to  the  local  boundary.  Hence  omitting  the  local  vorticity  (or  strength 
of  the  vortex  sheet)  at  the  point  where  the  B.C.  are  being  applied,  the 
remaining  components  of  velocity  must  sum  to  a  vector  parallel  to  the 
local  boundary.  Another  way  of  stating  this  condition  ’is  to  say  that 
if  all  components  of  the  velocity,  excepting  the  local  vorticity,  sum  to 
a  resultant  which  is  parallel  to  the  boundary,  the  subsequent  addition 
of  the  vorticity  component  will  not  affect  that  parallelism  (only  the 
magnitude) . 
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where  gg  is  the  so-called  single  airfoil  contribution 


% 


5 


[  A**  ^  (%)  dx-j  cvQ 


^  % 

F-* 


and  the  kernel  K  is  given  by 


A? 

r-* 


r?) 


(  Q 


If  a  modified  vorticity  function,  f,  is  defined  by 


V/W' 


Vi  L  - 


(/C) 


then  the  integral  equation  (7)  takes  the  alternative  form 


£(%) 


£(7))  K  (*;!])  Jyj 


( M ) 


In  order  to  obtain  numerical  solutions  to  this  equation  the  follow¬ 
ing  coordinate  transformation,  followed  by  discretization  of  the 
variables,  is  introduced. 
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/f\ 
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i 
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x  - 

k( 

)  -  cjrtG) 

,  Air* 

Lx 

f  l?#  1 
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\ 
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i 

i  - 

v-  iL 
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where  f/A  and  p  are  positive  integers  less  than  or  equal  to  n , 
or  else  zero. 


Under  this  parameterization  Eqs .  (8)  and  (9)  reduce  to 
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^te)-  4  6  (' 

K(W'  i  £b  -«3'H[-F(iatn4'-£  c«-i,V)  4-  4-)3  * 


.  itt 

O;  8 

and  the  integral  appearing  in  Eq.  (11)  reduces  to 


fe’j 


^  fo)  ^y)Ay 


r " 

j  <f  6^»v<l>  KM') 

'Jr) 


If  the  numerators  of  the  integrals  in  Eqs.  (8')  and  (9*)  are  assumed 
to  be  expanded  in  cosine  series,  it  may  be  shown  that  these  Glauert- 
type  integrals  may  be  approximated  by  sums 


/VU 

{a)*-*  Z a") 


rr^ 


where 


2+4) 

X  ttn 


m  ^tk 


±l*WH 


On- 


\-/vt 


and  S  is  the  Kronecker  delta, 
rs 


OirAfn 

f-A-o 

f  'J-A4 


04) 
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In  Eq.  (13),  on  Lhe  other  hand,  a  direct  application  of  the  trape¬ 
zoidal  rule  results  in 


i&i*  ( 4y0 


When  these  individual  results  are  combined  in  integral  equation 
(11)  there  results  a  set  of  n/1  simultaneous  equations,  each 
equation  corresponding  to  a  particular  value  of  the  index  J_  . 
These  may  be  expressed  in  matrix  notations  as 

A  ] 


Both  and|l/jare  column  matrices  with  n/1  rows  in  which  the  JL  ^ 
rows  are  given,  respectively,  by 


SU) 

9^  fit™  3 

>  -v  — ■  /"N 


/K 


(  \  ^ 
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Partitioning  u  into  two  parts 


U  —  o(  JJ, 


d-tr?  ol  Uj, 


where 


uT  =  *L*,t  i- 

7-0  7 


“TV 

„  /v  *  ‘ 

Ctrl  1  — 
/>V 


]  * 


rLxv\}~  ^ 

The  solution  of  Eq .  (15)  may  be  expressed 


r  i  tt) 


(-Jo) 
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J<  /F 
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(z  2) 


where 


,4"'  I/(° 


(  ^  3  , 


.  A"  U(* 


(2*) 


Cascade  Properties 
The  circulation  round  one  blade  is 

r  «  c  J  j (x)<k  -  cV\  -  c\ (1 

<5 


V  -  ic.VJ  (i A<une)  ^  (e>)de 
If"  ii(l)Z  o- 


n  ~ 

sV  ' 

X  - 
-\  /  “ 

Gl 

~"R.  t<rj  cC 


where 


/Tt 


Knowing  the  circulation  allows  the  mean  flow  velocity  V  and  direc¬ 
tion  o(  to  be  related  to  the  upstream  (or  downstream)  resultant 
flow.  Arranging  the  cascade  along  a  vertical  axis,  the  various 
velocity  vectors  take  on  the  aspect  shown  in  the  sketch. 


From  the  geometry  of  this 
sketch  we  deduce 


-tan  o(  =  tanod 
X  l 


2  tan  ) 


&  f  ' 


whereas  the  dimensionless 
circulation  for  one  pitch 
distance  along  the 
cascade  is  given  by 


_T  _ 
s  V 


£*?(•{+£)  []Wo/v~ 


Replacing 
Eq.  (25)  to  Eq. 


x  in 
(29) 


Eq.  (29) 
results 


by  Eq. 
in 


(28)  and  then  equating 


I  -f  ^  Cto  fo  —  AtAyi* 
4  Q  +  /Km  jl  4- 


15 


Presuming  thatch  may  be  determined  by  this  means,  it  is  then  possible 
to  express  V  in  terms  of  the  physically  determinate 


W,  Ori  °C> 

C-O-J  ( 


W,  c^cn  o(. 


f  /  i  <^r 


°<i 


//2-vh  f)j  Tam. 
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This  results  in  the  interesting  conclusion  that  a  given  cascade  may 
be  solved  up  to  the  determination  of  R  and  Q  without  regard  to  the  up¬ 
stream  flow  conditions.  Subsequently  the  velocity  triangles  may  be 
deduced  for  any  upstream  flow  condition  by  finding  \/ac  from  (30) 

and  (31)  and  applying  (28) .  In  fact  these  manipulations  may  be  used 
to  derive  Merchant's  formula 


""R  A  j >0 

2  fijit.  ft  H  Q  -Tt  tdM  yG 


-f 


X/Uc  T  -  G*.  +  R  h^/l  j 

-£— - - l&Y  ^ 

2  fu-c.  jrtJr  Q  -  ji 


showing 


the  linear  dependence  of 


(  y  2. 
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APPENDIX  B 


PLATE  MECHANICS 

The  moment  acting  on  the  spring  supported  plate,  see  Fig.  2,  is 
expressed  by  first  determining  the  aerodynamic  moment  about  the 
leading  edge  ,  r 

A10Mo  ~  <=■  fxctt  -  (&) 

/o  0 

This  allows  the  computation  of  the  rearward  spring  force 
immediately  from  equilibrium  considerations 

fr  -  AJ  /C*Vv-  (j&2) 

Then  computation  of  the  aerodynamic  lift 

lift  -  fotl  -  J 'fC'U')  (/  V')  (&3) 

allows  computation  of  the  leading  edge  spring  force 

■fC  =  lift  -  -fr 

Knowing  the  spring  constants  k£  and  kr  the  leading  edge  and 
rearward  support  deflections  may  be  computed  in  turn 

y(r,0 

yf  0  -  -fr/kr  (BS) 

where  Y  is  the  colum  matrix  (vector)  of  elastic  displacement 
dimensionless  w.r.t.  the  chord  and  the  indices  are  given  in  terms 
of  1  <  indix  <  m  s  n+1. 


The  bending  moment  distribution  may  now  be  computed  rearward 
from  the  trailing  edge  according  to  . 

At  —  fee-  X)  M  *  Vfa  -  c)  -fr  for  -X)  (yG) 


At~ ->-V'ySr-C)fr£xr-x) 


Co r&  —c-oj  £>r 

a 


where  U  is  the  unit  step  function  and  f ,  T  are  dummy  variables  of 
integration  analogous  to  x,  0. 


In  these  integrals  (Bl,  B6  and  others)  the  integrals  are 
determined  at  the  discrete  mesh  points;  hence  the  numerical 
integrations  may  be  evaluated  using  either  a  trapezoidal  rule,  a 
Simpson’s  rule  or  other  appropriate  numerical  technique. 


For  the  approximation  to  the  curvature  d2ydx2=cd2Y/dX2  several 
trial  models  were  evaluated.  Eventually,  in  order  to  perform  the 
integration  with  uniform  steps  in  the  angular  variable  (0  or  if)  the 
following  approximation  was  developed. 
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With  Y«y/c  and  X*x/c=h (l-cos  0}  it  is  found  that 
'll*-  =  ^  [sin.  9  ^  -  «“  <5> 

and  expressing  the  derivatives  in  terms  of  control  differences, 

c%>'3^ -J>  = 

\hs!^)Cy(e-i.O-2y(e,ij+y(e*hO-ihIcc>sCe^){-X(r-i,0*y(f'y))]/^^ 


where  u=»r/m  and  a  common  factor  h  is  retained  to  reduce  the  size  of 
the  elements  in  the  D  matrix. 

In  the  application  of  this  result  it  is  found  that  the  first 
and  last  rows  in  the  (n+l)x(n+l)  matrix  are  not  required  since  the 
bending  moment  is  zero  at  leading  and  trailing  edges.  Recognizing 
that  two  values  of  Y(0,1)  are  known  from  previous  equilibrium 
considerations,  two  columns  may  be  removed  and  the  useful  D  matrix 
may  be  reduced  to  (n-l)x(n-l)  dimensions.  In  this  reduced  form  the 
solution  for  the  remaining  Y  terms  may  be  obtained  by  inverting  the 
reduced  D  matrix  and  solving  for  the  remaining  unknown  values  of  Y 


where  the  subscript  r  refers  to  the  ^educed  versions  of  these 
matrices. 

Upon  reassembling  the  full  Y  vector  it  is  then  possible  to 
estimate  the  elastic  slopes  of  the  camberline  using  central 
differences*  to  correct  the  stagger  angle  and  the  camberline  slopes 
for  the  next  iteration 

^new=^oid+arctan  (*  (m ,  1 )  -  Y  ( 1 , 1 )  > 

Y*  ■  Z 

Z^^tanCarctan  Zeld  +  arctan  (Y  (m,  1) -Y  (1 , 1) ) 

Adding  the  elastic  slopes  to  the  unstrained  camberline  slopes  the 
total  slopes  are  then  made  available  as  inputs  to  the  aerodyamics 
subroutine  for  the  next  iteration. 


.For  leading  edge  and  trailing  edge  resp.  forward  and  backward 
differences  are  used. 
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PPENDIX  C 
Computer  Programs 
test  00. m 
solve. m 

in  MatLab  language 
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— 5 

Ye\-^trrr  a-*  e.  -tl 
tv  V*  U.'t’  **tU  v  Vx 


g.  3 


Computation  Flow  Chart  for  solve.ro  Program  (Mat lab) 


Rigid  and  Compliant  Blade  Comparison 


Rigid  Blade 


Compliant  Blade  -  onspeed 


Compliant  Blade  -  10%  overspeed 


Flow  Coefficient 


Leading-Edge  Loading  Parameter  for 
Rigid  and  Compliant-Bladed  Rotors 


Mi**  vp  • :  1  or»p  %%%%% 


Coefficients  for  rearward  force 


epdrArr/  (1  *(id*  l  pa  ram***e*  <  \  pa  ram-?  I  ♦  f  t  *  1 1  par  am"  31)  ; 
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